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■ ABSTRACT 

> 

O , 

Q , We consider the observational signature of the dynamical effects on the lumi- 

! nosity function of globular clusters. For the three best studied systems, in Milky 

CN ! Way, M31, and M87, there is a statistically significant difference between the in- 

ner and outer population of globular clusters. In all cases the inner clusters are on 
average brighter than the outer clusters (0.26 < Amo < 0.84) and have a smaller 
i dispersion in magnitudes (0.04 < Act < 0.53), with the larger differences for the 

local, better observed samples. The differences are of the type that would be 
expected if the inner population had been depleted by tidal shocks. The results 
^ I suggest that the inner population suffers substantial evolution from its initial 

^ I distribution and cannot not be used as a standard candle without correction for 

^1 dynamical evolution. 

O 

■ Subject headings: globular clusters: general — galaxies: star clusters — galaxies: 



individual (M31, M87) 



1. Introduction 

The turnover magnitude of the luminosity function of globular clusters in external galax- 
ies has been used as a standard candle for distance measurements (cf [Jacoby et al. 199^ ). 
The method is usually applied to elliptical galaxies with the calibration to the globular clus- 
ter system (GCS) in the Local Group. Although it seems to be quite universal, the center 
of the luminosity function can be expected to differ from galaxy to galaxy. In fact, there 
is no consistent theoretical model that predicts the same mean mass or luminosity for the 
initial distribution of globular clusters (but see [Fall fc Rees 1985| , and [Vietri fc Pesce 19"95| ). 



Recent work by [Gnedin fc Ostriker (1997")| and [Murali fc WeinbergI (1996a,b), which built on 



less detailed analysis of [Aguilar, Hut, fc Ostriker (1988)| and |Ghernoff fc Weinberg (1990) 



demonstrated that the GCS in our Galaxy suffers a substantial depletion over the Hubble 
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time due to a variety of dynamical processes, but primarily the tidal bulge and disk shocks 
and the enhanced evaporation of stars through the tidal density cutoff. These effects vary 
with distance to the Galactic center and are most pronounced in the inner part. In this 
Letter we look for an observational evidence for the imprints of the dynamical effects and 
address the question of universality of the globular cluster luminosity function (GCLF). 

If the dynamical effects are important, the inner population should show a distribution 
different from that of the outer population that presumably has the initial form (modulo 
internal two-body relaxation process). For our Galaxy three-dimensional information is 
available, but for others only projected radii are known. Therefore, the simplest approach 
is to divide the observed sample into two halves, inner and outer, according to the distance 
from the center of the galaxy. 

The GCLF is conventionally fitted to a Gaussian function in magnitude, 0(m) oc 
exp (— (m — mo)^/2(T^). It provides a useful measure of the overall shape and the parame- 
ters (mo, a) of the GCS. We assume also that the two halves of the sample can be fit by 
the same function, but with different parameters. This allows us an easy way to compare 
the two populations. As an independent check that does not use the assumption of normal 
distribution, we compare the median points (/i) of the two populations. 

We use the maximum likelihood estimator (MLE) as a prime tool to evaluate the param- 
eters of the distribution and their uncertainties ( [Lupton 1993|) . |Hanes fc Whittaker (1987] 



reviewed this method in detail and formulated the necessary corrections for the magnitude- 
limited data. We use their equations (lOa-lOc) to estimate the mean and dispersion of the 
sample as well as the normalization coefficient, all of which are allowed to vary. [Seeker (1992) 



investigated a number of functional forms to be applied in the ML analysis of the GCLF 
and found the Student's t-distribution to be optimal. Since the deviation in the inferred 
parameters as compared to the Gaussian is small, we chose the latter as it is simpler and 
more commonly used. In most studies up to now the turnover point and the dispersion of the 
GCLF were determined by fitting a Gaussian to the histogram of the magnitude distribution. 
This method implies binning of the data and leads to a significant loss of information that 
becomes especially important in case of a small sample (as in the Milky Way and M31). We 
use it as a secondary estimator, mainly for consistency with the previous results. 

Monte-Carlo simulations were performed to establish the statistical significance of the 
differences of the inner and outer populations. The original sample was divided into two 
parts randomly, and for each such realization the differences of the means, Amg, and of 
the medians, A/i, were recorded. 10,000 of the random realizations constitute a very nearly 
gaussian sample from which the probability of drawing a given Amo or A/i is calculated. 
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2. Milky Way 

The sample of the Galactic globular clusters is best studied and essentially complete. We 
take the recent compilation by pjorgovski (1993"}] of 140 clusters with measured magnitudes 
and distances from the Galactic center. Sorted by the distance sample is divided into two 
equal parts at the boundary radius of R1/2 = 6.6 kpc. Note that |Zinn (1993)| used a similar 
division to separate the old and young halo clusters. 

The estimated parameters of the Milky Way GCLF are given in Table ^ The mean 
of the whole distribution is at My = —7.17 ± 0.12, in accord with Seeker's (1992) estimate 
of = -7.14 (fitting the histo gram, however, gives —7.36 ± 0.17; [Harris et al. 1991| ). 



The means of the two populations differ by almost 0.7 mag, and the medians by 0.33 mag, 
with the inner part being systematically brighter. The outer clusters also have much larger 
dispersion. The result is statistically significant as given by the Monte-Carlo probability 
of less than 0.5% (Table |l|) that the observed ArriQ occurs by chance. As expected, the 
dispersion of Amo from the Monte-Carlo realizations is nearly the same as the MLE of the 
standard error of Amo. The difference in medians is smaller but nonetheless noticeable. We 



have checked that restriction of the sample to the halo clusters (with [Fe/H] < —0.8; [Zinn 



19851) changes all parameters by much less than their standard errors. 

The difference of the inner and outer populations is illustrated on Figure |l] that shows 
histograms for the two subsamples and our fit obtained with the ML estimated parameters. 
Fitting a Gaussian directly to the histograms depends on the size of the bins. We found a 
slight variation within the errors of the two means obtained this way for a range of bin sizes, 
dm, from 0.3 to 0.5 mag. Table ^ shows the results for dm = 0.45 mag which are consistent 
with the ML estimates. The dispersion of the outer sample is less well constrained; note also 
the large standard error of its ML estimate. The difference of the means is reduced to 0.5 
mag, but still is a two-sigma result. 

Kavelaars & Hanes (1996)' have recently undertaken a similar study and found no dif- 
ference in the centers of the inner and outer populations of the halo clusters, though the 
dispersions were significantly different (and close to our estimates). To facilitate a compar- 
ison with their results, we tried to match the sample they used as close as possible; it was 
drawn from the database maintained by W. HarrisQ. The second half of Table |1] shows that 
Amo is still significant at 1.5-sigma level when calculated using ML, but disappears when 
inferred from the histogram fitting. For the latter we used the same bin size (0.4 mag) as 



Kavelaars fc Hanes (1996)| . We doubt the trustiness of the latter estimate based on binning 
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Fig. 1. — Histograms for the inner (solid) and outer (shaded) populations of the globular 
clusters in Milky Way. The two Gaussians with the ML estimated parameters are superim- 
posed on the histograms for comparison. The solid line is for the inner sample, the dashed 
is for the outer. The peak of the inner population is brighter by 0.7 mag, while the outer 
population has considerably larger dispersion. 
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of already very small sample (A^ < 50). The lesser statistical significance of the ML estimator 
of ArriQ in this case arises probably from the smaller size of the sample. 



3. M31 

The Andromeda has more globular clusters than Milky Way, and the dust obscuration 
is less of a problem. We use the most recent compilation of the GCS in M31 by J.-M. 
PerelmuterQ with 1 cluster rejected on basis of its radial velocity (the kinematic data were 
kindly provided by J. Huchra). The sample becomes incomplete at V ^ 17.5, well beyond 
the peak of the distribution. We considered, therefore, two cuts of the original data at 
^ = 18 and V = 17.5, respectively. 

Figure |^ shows the histograms for the more conservative second cut. The centers of 
the two populations are clearly different with Amo ~ 0.8. The difference in dispersions is 
smaller than in Milky Way and is within the errors. 

The ML estimates agree very well with the histogram fitting results (see Table |^) for both 
cuts of the sample. The inner population is consistently brighter than the outer with huge 
statistical significance (7-sigma for Amo and 6-sigma for Ayu). Thus there is little doubt that 
dynamical evolution left a strong imprint in the inner parts of the galaxy. Relative closeness 
of M31 makes the completeness limit much better off than in M87, so the result is essentially 
unaffected by the observational selection. 



4. M87 

The giant cD galaxy at the center of the Virgo cluster possesses the largest sample of 
globular clusters observed around galaxies, although the completeness is a severe problem. 
The most recent data available to us is from [McLaughlin, Harris, fc Hanes (1994)| . All 



sources have V < 24, and the completeness limit is about V ~ 23.5. Along with the real 
clusters there could be a contamination from the foreground stars or background sources. 
We restricted the sample to a range of radii (1.'21 < R < 7') where the signal-to-noise ratio 
is greater than unity. 

Figure ^ shows the histograms for the sample cut at V = 24. The inner and outer 
populations have similar distributions but there is an apparent excess of the outer clusters 
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Fig. 2. — Same as Figure |I|, but for the globular clusters in M31. The sample is restricted 
to V < 17.5 due to incompleteness of the data. 
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at the faint end. Also, the inner clusters prevail for 21 < < 23 (but not in the very bright 
end). 

The MLE of the difference of the means is Amo = 0.26 ± 0.03 (Table D and is formally 
very significant. These standard errors should be taken with caution, however, as the sys- 
tematic errors due to incompleteness of the data are likely to be larger. For example, the 
difference of the medians is only 0.16 mag. Similar results are obtained when the sample is 
cut a.tV = 23.6, i.e. very close to the peak of the GCLF. In both cases statistical corrections 
applied to ML estimates of mo and a are important. Therefore the results are somewhat 
sensitive to the completeness limits. A deeper sample is required for robust comparison of 
the inner and outer parts of the GCS. 

Fitting the histograms becomes unstable because of the incompleteness. Thus, even 
for the whole sample [McLaughlin, Harris, fc Hanes (1994)| found that a broad range of mo 



and a gives similar "goodness" (x^) of the fit. Our fits give fainter central magnitudes and 
broader dispersions, but the scatter is very large. Also, different parameters were obtained 
for different bin sizes. For this reason we consider the fits as unreliable. 



5. Conclusions 

We investigated the differences of the inner and outer populations of globular clusters 
in Milky Way, M31, and M87. There is a considerable difference in the turnover magnitudes, 
Amo ~ 0.7 — 0.8, for the two close and more complete samples. Monte-Carlo simulations 
confirmed statistical significance of the result. An independent check is provided by using 
the medians instead of the means which is unaffected by the wings of the distribution. 

In all three galaxies, the inner cluster population has brighter mean magnitude and 
smaller dispersion. Presumably, the differences arise due to dynamical processes operating 
most efficiently in the inner part of the galaxies. The sign of the observed effects is as 
expected for a depletion in the inner parts of the galaxies of low luminosity, low density 
clusters. We will investigate the dynamical interpretation of the results in a complementary 
study ( pstriker fc Gnedin 19961 ). 

If the dynamical changes in the GCLF are significant, it cannot be used directly as 
a standard candle. Special corrections should be introduced that take into account the 
evolution. 



I am greatly indebted to Prof. J. P. Ostriker who proposed this investigation and 
provided constant support and encouragement. I would like to thank R. Lupton for the 
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Fig. 3. — Same as Figure |I|, but for the globular clusters in M87. Only clusters with V < 24 
are shown, and the completeness limit is about V = 23.5. 
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numerous help with the statistics. I am grateful to W. Harris, J. Huchra, D. McLaughlin, 
and J.-M. Perelmuter for providing the data on the clusters in M31 and M87. This project 
was supported in part by NSF grant AST 94-24416. 
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Table 1. Globular Clusters in Milky Way 



Sample 


N 


mo (mag) 


a (mag) 


/i (mag) 


All clusters 


Total 


140 


-7.17±0.12 


1.43 ±0.49 


-7.29 ±0.09 


Inner 


70 


-7.51 ±0.13 


1.08 ±0.44 


-7.40 ±0.12 


Outer 


70 


-6.82 ±0.20 


1.65 ±0.68 


-7.07 ±0.22 


difference 




0.69 ±0.24 




0.33 ±0.26 


probability*^ 




4.7 X 10^3 




9.7 X 10-2 


Fit (bins 0.45 mag) 


Inner 




-7.48 ±0.16 


1.04 ±0.07 




Outer 




-6.99 ±0.24 


1.44 ±0.12 





Comparison with [Kavelaars fc Hanes (1996' 



Total 

Inner 

Outer 

difference 

probability*^ 



93 
49 
44 



-7.16 ±0.15 
-7.36 ±0.15 
-6.94 ±0.26 
0.42 ± 0.30 
0.15 



1.40 ±0.53 
1.02 ±0.46 
1.71 ±0.79 



-7.35 ±0.11 
-7.43 ±0.07 
-7.25 ±0.32 
0.18 ±0.33 
0.33 



Fit (bins 0.4 mag) 


Inner 


-7.38 ±0.12 


0.64 ±0.10 


Outer 


-7.39 ±0.32 


1.39 ±0.20 



^^Probability of the Monte-Carlo realizations that a given result occurs by chance. 

'^Sample constructed to match that used by [Kavelaars fc Hanes (1996]] from the database 
of W. Harris. The clusters have [Fe/H] < —0.8 and lie within 1 < i? < 85 kpc from the 
Galactic center. The dividing boundary is at R1/2 = 8 kpc. 
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Table 2. Globular clusters in M31 



Sample 


N 


mo (mag) 


a (mag) 


1^ (mag) 


Clusters with V < 18 


Total 


168 


16.29 ±0.06 


0.88 ±0.27 


16.32 ±0.10 


Inner 


84 


15.93 ±0.08 


0.70 ±0.27 


16.01 ±0.07 


Outer 


84 


16.71 ±0.09 


0.94 ± 0.32 


16.68 ±0.09 


difference 




0.78 ±0.12 




0.67 ±0.11 


probability^ 




4.1 X 10-1° 




7.3 X 10-^ 


Fit (bins 0.3 mag) 


Inner 




15.95 ±0.10 


0.73 ±0.14 




Outer 




16.70 ±0.12 


0.78 ±0.17 




Clusters with V < 17.5 


Total 


160 


16.32 ±0.06 


0.88 ±0.25 


16.25 ±0.12 


Inner 


80 


15.93 ±0.08 


0.71 ±0.27 


16.01 ±0.07 


Outer 


80 


16.76 ± 0.08 


0.95 ± 0.30 


16.63 ± 0.08 


difference 




0.84 ±0.11 




0.62 ±0.11 


probability^ 




4.7 X 10-12 




1.9 X 10-^ 


Fit (bins 0.35 mag) 


Inner 




15.88 ± 0.09 


0.67 ±0.16 




Outer 




16.83 ± 0.26 


0.93 ±0.13 





^Probabihty of the Monte-Carlo reahzations that a given result occurs by chance. 
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Tabic 3. Globular Clusters in M87 



Sample 


N 


mo (mag) 


a (mag) 


1^ (mag) 


Clusters with V < 2A 


Total 


3287 


23.22 ±0.02 


1.26 ±0.15 


22.92 ±0.02 


Inner 


1644 


23.09 ±0.02 


1.24 ±0.18 


22.84 ±0.03 


Outer 


1643 


23.35 ±0.02 


1.28 ±0.18 


23.01 ±0.03 


difference 




0.26 ±0.03 




0.16 ±0.05 


probability*^ 




4.1 X 10"^^ 




2.9 X 10"^ 


Clusters with V < 23.6 


Total 


2718 


23.06 ±0.02 


1.21 ±0.15 


22.69 ±0.02 


Inner 


1359 


22.93 ± 0.02 


1.18 ±0.18 


22.63 ± 0.03 


Outer 


1359 


23.19 ±0.03 


1.24 ±0.18 


22.73 ± 0.03 


difference 




0.26 ± 0.04 




0.10 ±0.04 


probability^ 




2.4 X 10-^3 




1.1 X 10-2 



^Probabihty of the Monte-Carlo reahzations that a given result occurs by chance. 



